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Abstract. Low-Mass X-ray Binaries (LMXRBs), believed to be tiie progenitors of recycled millisecond pulsars (MSPs), 
occur abundantly in globular clusters (GCs). GCs are therefore expected to host large numbers of MSPs. This is also confirmed 
observationally. The MSPs continuously inject relativistic electrons into the ambient region beyond their light cylinders, and 
these relativistic particles produce unpulsed radiation via the synchrotron and inverse Compton (IC) processes. It is thus 
possible, in the context of General Relativistic (GR) frame-dragging MSP models, to predict unpulsed very high energy 
radiation expected from nearby GCs. We use a period-derivative cleaned sample of MSPs in 47 Tucanae, where the effects of 
the cluster potential on the individual period derivatives have been removed. This MSP population is likely to have significant 
pair production inhibition, so that slot gaps and outer gaps are not expected to form in the pulsar magnetospheres. The 
utilisation of unscreened pulsar potentials is therefore justified, and fundamental tests for the predicted average single pulsar 
gamma-ray luminosities and associated particle acceleration are simplified. Using a Monte Carlo process to include effects 
of pulsar geometry, we obtain average injection spectra (with relatively small errors) of particles leaving the MSPs. These 
spectra are next used to predict cumulative synchrotron and IC spectra expected from 47 Tucanae, which is a lower limit, as 
no reacceleration is assumed. We find that the IC radiation from 47 Tucanae may be visible for H.E.S.S., depending on the 
nebular field B as well as the number of MSPs A' in the GC. Telescopes such as Chandra and Hubble may find it difficult to 
test the SR component prediction of diffuse radiation if there are many unresolved sources in the field of view. These results 
may be rescaled for other GCs where less information is available, assuming universal GC MSP characteristics. 
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INTRODUCTION 

A total of 137 globular cluster (GC) pulsars have been 
discovered in 25 GCfl following the discovery of the 
first GC milUsecond pulsar (MSP) in M28 tl . Low-Mass 
X-ray Binaries (LMXRBs), believed to be the progeni- 
tors of recycled MSPs [2], occur abundantly in GCs. GCs 
are therefore expected to host large numbers of MSPs, up 
to ~ 200 MSPs or more [|3i, i]. Indeed, GC MSP spin 
properties seem consistent with the recycling scenario 

Terzan 5, 47 Tucanae, and M28 collectively contain 
nearly half of all GC pulsars, housing 33, 23, and 1 1 pul- 
sars respectively f?]. GC MSPs are sources of relativistic 
electrons, which are continuously being injected into the 
ambient region beyond the MSPs' light cylinders. These 
relativistic particles produce high-energy emission via 
synchrotron radiation (SR) and inverse Compton scatter- 
ing (ICS) on bright starlight photons as well as on the 



http://www.naic.edu/~pfreireGCpsr.html 



cosmic microwave background (CMB). 

In this paper, we calculate the cumulative injection 
spectrum and resulting unpulsed SR and ICS fluxes, us- 
ing a population of 13 MSPs in 47 Tucanae, with cor- 
rected values of their period-derivatives GC gamma- 
ray visibility is also discussed. We use a more refined 
injection spectrum, originating from a General Relativis- 
tic (GR) frame-dragging MSP model (e.g. |jlS[i]), than 
that assumed by |4], and our calculations are comple- 
mentary to pulsed gamma-ray flux predicted by ifToll . As 
we only consider particles originating from MSP magne- 
tospheres, with no further acceleration, our calculations 
should be viewed as lower limits to the expected TeV 
flux. 



INJECTION SPECTRUM 
CALCULATION 

As in lH^lTl, we use a population of 13 MSPs in 47 Tu- 
canae, with corrected values of their period-derivatives 
P J^. We calculate the injection spectrum Q' of elec- 
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FIGURE 1. Loss timescales (left panel for B = 1 flG, right panel for 6 = 10 flG). Solid lines represent Z = 0, dashed lines Z = 1 , 
thin lines Tj^ij, intermediate straight lines Tesc, and thick lines Teff (see text for details). 



Irons leaving each MSP /, with / = I,- -,13, by bin- 
ning the number of primary electrons leaving a stellar 
surface patch and moving along a B-line according to 
E^'^ = yic'iieC^, the residual electron energy at the light 
cylinder, divided by energy bin size. 

We next randomly choose = 100 MSPs (with ran- 
dom incUnation angles x), and sum their particle spectra 
to obtain a million cumulative spectra from the Monte 
Carlo process: 



6; 



cum 



k=l 



dEf, 



(1) 



where for each index j = 1 , 2 , . . . , M = 1 0^, a total of N = 
100 particle spectra (randomly sampled from 104 spec- 
tra: 13 MSPs X 8 values of x, X ^ 10°, 20°, • • • 80°) are 
summed. We therefore oversample from 13 to 100 pul- 
sar particle spectra, to obtain the cumulative injection 
spectrum for 100 MSPs into the GC 47 Tucanae. 
Note that the single MSP particle spectra Q' are not func- 
tions of observer angle ^, so that the relative errors of 
the average cumulative particle spectrum < Qium > are 
smaller than for the average cumulative gamma-ray spec- 
trum < {dNy/dE)iam > at earth, which was calculated in 
flQl . The relative uncertainty of < Qcum > is due to the 
different values of P and P, as well as different inclina- 
tion angles x, of the MSP population members. The thick 
solid lines in Figure|2]indicate < Qium >■ 



UNPULSED GAMMA-RAY FLUX 

We divide the region where unpulsed radiation is gener- 
ated into two zones: 'zone 0' (Z = 0) reaching from r — Q 
to r = rc, with the core radius, and 'zone V (Z — 1), 



reaching from r = rc to r = r^^, with rhm the half mass 
radius. Using the cluster core formula (Eq. 5 in we 
find an energy density Wrad 3000 eV/cm^ for Z = 0, 
and Mrad ~ Loc/i'^TtT^c) ^ 100 eV/cnr' for Z = 1 for 
the bright starlight component (corresponding to temper- 
ature T « 4500 K), with Lqc the observed cluster lumi- 
nosity, and 7 = (rc + rhm) /2. We use similar parameters 
as |4]: distance d = 4.5 kpc, rc = 0.58 pc, r^m — 3.65 pc, 
and Lgc = 7.5 x lO^L©. For the CMB component, we 
use McMB ~ 0.27 eV/cm^ (corresponding to T « 2.76 K). 

The radiation loss, escape, and effective timescales 
(Figure [Til are given by 
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The left panel of Figure [T] is for a nebular field B ~ 
1 ;UG, while the right panel is for B = lO^uG. Solid Unes 
represent Z = 0, and dashed lines Z = 1 . 

For B = lOjtiG, Tesc B is 10 times larger than for 
B =1 /iG. The bright starlight component dominates Trad 
at low energies for all B-field strengths, while the CMB 
component dominates for B — I fiG at high energies, 
and Tiad is ^ 10 times lower for the case of B = 10 jJ-G 
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Bottom thick lines: injection spectra Q x 10^ (/TeV/s); thin hnes: Teff x 10^" (s); top intermediate hnes: steady-state 



FIGURE 2 

electron spectra dNf/dE,, (/TeV) (left panel for B 
Z — I. We used dNg/dEg =< >cum x feff 



1 /xG, right panel for B = lO^G). Solid lines represent Z = 0, dashed lines 



than for the case of Z? = 1 /iG at high energies. When 
B = I /xG and B ~ 10 ^G, Trad ^ fesc at small Eg, while 
Trad ^ '^esc at large Eg for B — IjxG. Thus Tgff ~ T^d at 
small Ef. for all field strengths, while Tgff « Tesc at large 
Ee for B = l /xG, but Tgff « Trad for B = 10 ^G and Z = 1 . 
The steady-state electron spectrum is calculated using 



-< 2cum > 't'eff 



(7) 



for each zone. Figure |2] indicates the (scaled) injection 
spectrum for Z = by thick solid lines (average of Eq. 1). 
We divide dN/ dE^ of Z = by Tesc of Z = 0, and obtain 
the thick dashed lines for the injection spectrum of Z = 
l.The product of Tgff (thin lines indicating scaled Tgff) 
with the injection spectra (thick lines) give the steade- 
state electron spectra for each zone (lines of intermediate 
thickness). As previously, the left panel of Figure|2]is for 
B = \ liG, while the right panel is for B = 10/xG. SoHd 
lines represent Z = 0, and dashed lines Z—l. 



DISCUSSION AND CONCLUSIONS 

The predicted SR and ICS fluxes are larger for B = 10/xG 
than for B = 1 jiG (the SR fluxes differ by about 2 orders 
of magnitude). The SR component occurs in the optical 
/ X-ray waveband, and the ICS (scattering starlight and 
CMB), in the gamma-ray waveband. 

Our predictions for the IC flux are notably smaller than 
those of [4]. This is mainly attributed to the fact that 
^ assumed an average spindown luminosity per MSP 
of ^ 10^^ erg/s, while we find ^ lO^'* erg/s when us- 
ing the population of 13 MSPs with known 'cleaned' 
P (although we both assumed = 100 GC members). 



Also, pJ] used an average particle conversion efficiency 
of rj — 0.01, which is quite close to our numerical re- 
sult of T] = 0.0074 obtained in |10]. Lastly, our spectral 
shapes differ, as |4j] approximated the injection spectra 
by power laws with different spectral indices and cut- 
offs, while we calculated the cumulative injection spec- 
trum using the GR frame-dragging model as well as an 
actual GC population of MSPs. 

Due to the very bright starlight component in the 
GC core, the ICS from 47 Tucanae may be visible for 
H.E.S.S., depending on the nebular field B as well as the 
number of MSPs in the GC. (A lowest limit may be 
obtained by multiplying the graphs by 0.23, thereby as- 
suming = 23). Telescopes such as Chandra and Hub- 
ble may find it difficult to test the SR component pre- 
diction of diffuse radiation if there are many unresolved 
sources in the field of view. 

The results above 1 TeV may be roughly scaled for 
other GCs, where less information is available, by the 
factor 

100 y ^ d J ^(mtuc) 



(8) 



with < M >^ LGc/(4?rr^rn)' assuming universal GC 
MSP characteristics. As only MSPs were considered as 
sources of relativistic particles, and no further particle ac- 
celeration has been assumed, our estimates for the TeV 
flux should be viewed as lower limits. 
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FIGURE 3. Unpulsed SR (thin lines) and ICS (tliick lines) components (solid lines: B = I flG, dashed lines: 6 = 10 p.G; 
summation of Z = and Z = 1 spectra), with the latter consisting of a bright starlight and CMB component. Also shown are a 
prediction from |^ (short-dashed line), and the H.E.S.S. sensitivity 1 13]. 
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